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A sensitive search for the rare decays —> Att” and H® ^ pTr~ has been performed using data 

from the 1997 rnn of the HyperCP (Fermilab E871) experiment. Limits on other such processes do 
not exclude the possibility of observable rates for |AS'| = 2 nonleptonic hyperon decays, provided 
the decays occur through parity-odd operators. We obtain the branching-fraction limits B{Ll~ —>■ 

Att”) < 2.9 X 10“® and B(H® pTv~) < 8.2 x 10“®, both at 90% confidence level. 


The standard model allows IAS”! = 2 transitions 
through second-order weak interactions. This approach 
successfully describes AT® AT® mixing, which is currently 
the only observed IAS”! = 2 transition. While other 
AS'I = 2 processes have generally been considered too 
highly suppressed to be observed experimentally, it has 
been noted that the rate of AT® A'® mixing does not ex¬ 
clude nonleptonic |AS| = 2 hyperon decays at observable 
rates, provided that they proceed through new parity- 
odd channels [1, 2]. Measurements can thus be used to 
set limits on parity-odd contributions to hyperon decays. 
It is of interest, therefore, to perform sensitive searches 
for such decays. There is also interest in searches for 
direct jAS”! = 2 transitions in R-meson decays [3]. 

Observation of IAS”! = 2 nonleptonic hyperon decays 
at current levels of sensitivity would strongly suggest 
new physics. He and Valencia [2] have parametrized the 
strength of any new parity-odd interaction as a ratio to 
the strength of the electroweak interaction. This ratio is 
constrained by hyperon branching ratios [2], e.g. , 

H(S®^p7r-) = 0.9f^V : (1) 

\aEW/ 

where Onew (c^Ew) is the strength of the new (elec¬ 
troweak) interaction. Gurrent experimental limits on 
|AS| = 2 decays [4] include B{n~ —> Att”) < 1.9 x 
10“'® [5] and H(S® —> p7r“) < 3.6 x 10“® [6], both at 
90% confidence level (G.L.). We report a search for these 
decays using data from the 1997 run of HyperCP (Fermi¬ 
lab Experiment 871). These data are well suited for such 


studies as they contain large numbers of charged hyperon 
decays, ^10® S“ ^ Att” and ~ 10'’ fl~ —> AK~ decays. 

The HyperCP spectrometer is shown schematically 
in Fig. 1. The spectrometer (described in detail else¬ 
where [7]) was designed to have large acceptance for the 
decay chain S“ —> Att~ , A —> pTT~. In brief, a nega¬ 
tively charged secondary beam was formed by the inter¬ 
action of 800 GeV/c primary protons from the Tevatron 
ina0.2x0.2x 6 cm® copper target, with the secondaries 
sign- and momentum-selected by means of a 6.096-m- 
long curved collimator within a 1.667 T dipole magnetic 
field (Hyperon Magnet). The mean secondary momen¬ 
tum was about 160 GeV/c with momenta ranging from 
120 to 250GeV/c. Hyperon decays occurring within a 
13-m-long evacuated decay pipe were reconstructed in 
three dimensions using a series of high-rate, narrow-gap 
multiwire proportional chambers arrayed on either side 
of a pair of dipole magnets (Analyzing Magnets). 

The trigger for the data acquisition system [8] used 
scintillation-counter hodoscopes located sufficiently far 
downstream of the Analyzing Magnets so that the hy¬ 
peron decay products were well separated from the sec¬ 
ondary beam. A coincidence was required of at least one 
(“same-sign”) hodoscope hit consistent with a charged 
particle of the same sign as the secondary beam and at 
least one “opposite-sign” hit. To suppress muon and low- 
energy backgrounds the trigger also required a minimum 
energy deposit in the Hadronic Galorimeter. 

We searched for events consistent with either AS = 2 
decay chain, —> Att”, A —> pTT~ or —> S®7r“, 

S® —> pTT~. Also studied were the copious events from 
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FIG. 1: Plan view of the HyperCP spectrometer. Note that 
the z scale is compressed by a factor of 10 with respect to the 
X scale. 



FIG. 2: Event topology for all decays considered here. 




FIG. 3: The pn~n~ invariant-mass distribution prior to event 
selection (top) and from a Monte Garlo simulation of the 
fl~ —> Att”, a —> pTT~ process (bottom). 


the AK~, A ^ pTr~ decay sequence used for nor¬ 

malization. Such events all have the topology shown 
in Fig. 2, with three charged tracks forming two sepa¬ 
rated vertices. A least-squares geometric fit determined 
the positions of the primary and secondary vertices, as 
well as the chi-square (y^) value for the event to have 
the required topology. The reconstructed 12“ trajectory 
was traced back through the Hyperon Magnet using the 
measured magnetic field to determine its x and y coordi¬ 
nates at the midpoint of the target. Figure 3 compares 
the p'K~TT~ invariant-mass distribution of all events be¬ 
fore event selection with that of Monte Carlo-generated 
12“ —> A7r“, A —> pK~ events. 

Selection criteria for signal events were based on Monte 
Carlo simulations and studies of a purified sample of 
12“ —> AK~, A —> pTT~ events. Events were required 
to have three charged tracks with two tracks on the 
“same-sign” side of the spectrometer and one on the 
“opposite-sign” side. Most K~ events were 

excluded by requiring the invariant mass (treating all 
three charged tracks as pions) to exceed 0.500 GeV/c^, 
3 standard deviations (cr) above the K~ mass [4]. Since 
many background events originated from secondary in¬ 
teractions near the exit of the collimator or at the exit 
windows of the Vacuum Decay Region, all events were re¬ 
quired to have a primary-vertex z position between 150 
and 1180 cm from the exit of the collimator, well within 
the Vacuum Decay Region. The secondary vertex was 
required to be downstream of the primary vertex. The 


X^/degree-of-freedom from the geometric fit was required 
to be less than 1.7; this requirement was 98% efficient for 
12“ ^ AK~, A —> pTT~ events. Selection criteria for the 
projected 12“ position at the target were determined us¬ 
ing the purified 12“ —> AK~, A —> pTT~ sample. Due 
to differing resolutions in x and y at the target, an el¬ 
liptical target cut was used which was 92% efficient for 
12“ —> AK~, A —> pTT~. 

Additional selection criteria were specific to each 
mode. For the 12“ —> A7r“, A ^ p'K~ search we re¬ 
quired the invariant mass of the pn combination forming 
the secondary vertex to be within ±2.0MeV/c^ (±1.5 
times the rms mass resolution) of the A mass [4]. Decay 
polar angles in the respective hyperon rest frames were 
calculated with respect to the parent’s lab-frame momen¬ 
tum vector. Background events tended to have a small 
decay polar angle for the pion in the A7r“ center-of-mass 
system (0^ < 0.82rad), thus we required > 0.82rad. 
The invariant-mass distribution, assuming a ptt~tt~ fi¬ 
nal state, after all selection cuts is shown in Fig. 4a. No 
events were observed within ±9 a of the expected mass. 

For the 12“ ^ S*^7r“, —> p'K~ search, we required 

that the proton momentum be more than 38% of the 
p'K~'K~ total momentum, that the invariant mass of the 
pTT combination forming the secondary vertex be within 
±5.0MeV/c^ (±3.0 cr) of the mass [4], that the z po¬ 
sition of the daughter-hyperon decay vertex be within 
the Vacuum Decay Region, and that the polar angle of 
the proton in the pn center-of-mass frame be less than 
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FIG. 4: The pn~Tv~ invariant-mass distributions for events 
satisfying all selection criteria for (a) 14“ —> Att”, A —> p-K~ 
and (b) Q~ —> S°7r“, p'K~. 


TABLE I: Relative spectrometer acceptances vs. a decay 
parameters. 
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0.0 
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0.25 
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0.09 

1.0 
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0.50 

0.642 

0.87 

0.23 

- 1.0 

0.92 
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0.642 

0.79 

0.23 

0.0 

0.95 

1.00 

0.642 

0.70 

0.23 
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2.97rad. The invariant-mass distribution for 

events satisfying all selection criteria is shown in Fig. 4b. 
Again, no events were observed within ±9cr of the ex¬ 
pected mass. 

The spectrometer acceptances for these decays were 
estimated using a Monte Carlo simulation. The gener¬ 
ated momentum and position distributions at the 
target were tuned to match those observed in the data 
for a purified sample of Sl“ ^ AK~, A —> pTr~ decays. 
For all decay modes, the acceptance includes the prob¬ 
ability that the parent decay occur within the Vacuum 
Decay Region. The geometric acceptance of the spec¬ 
trometer was determined primarily by the apertures of 
the Analyzing Magnets. Due to the larger Q values 
for ^ Att” and ^ pTT~ compared to those of 
n~ —» AK~ and A —> p7r“, tracks from signal-mode 
decays would be approximately three times more likely 
to miss the magnet apertures and be lost. For events 
accepted in the magnet apertures, the trigger efficiency 
ranged from 93% for ^ Att", A —> pTT~ candi¬ 
dates to 99% for ^ AK~, A pTr~ events. The 
candidate signal- and normalizing-mode events were all 
accepted by the same trigger, so it was not necessary to 
cross-calibrate trigger efficiencies by mode. Offline event- 
selection efficiencies differed considerably by mode, since 
restrictive selection criteria to suppress background were 
required for the signal modes but not for the normalizing 
mode. To study the systematic uncertainties of the selec¬ 
tion efficiencies, the selection criteria were varied. This 
produced only slight changes in the results. 

The a decay parameter for ^ AK~ has recently 
been precisely measured and found to be small, but 
most likely not zero [9], while the measured value for 
^ S°7r“ is less precise and is consistent with zero [4]. 
Theory predicts little parity violation in the dominant 
decays [10] (as observed); however, we cannot ex¬ 
clude large decay asymmetries for —> Att~ . Sim¬ 
ilarly, we do not know the size of parity violation in 
^ pTT~. We therefore assign zero to the value of 
a for ^ Att” and S'’ ^ p7r~. The dependence of the 


acceptance on a is tabulated in Table I. The acceptance 
times the selection efficiency (including the maximal ac¬ 
ceptance variation due to uncertainty in a) was 
for > S^TT”, S® —> pn~, 6.9^2 2 ^ for ^ Att", 

A ^ _p7r“, and (34.9 ± 1.2)% for ^ AK~, A pTr~. 

Signal-mode branching ratios were normalized using 
the observed ^ AK~, A —> p7r“ sample, whose 
invariant-mass distribution is shown in Fig. 5. The fig¬ 
ure includes a fit of the mass distribution to a Gaus¬ 
sian plus a second-order polynomial. The number of 
normalizing-mode decays was estimated using three dif¬ 
ferent background-subtraction methods. In all three, an 
estimate of the background was subtracted from the to¬ 
tal number of events observed within ±3 a of the D mass. 
The first background estimate was the sum of all events 
within ±3 a of the D mass minus the integral of the fitted 
second-order polynomial over that region. In the second 
method, a fit was performed over a limited region (±3 a) 
around the D mass to a Gaussian plus a constant. The 
third method averaged the two bins at ±3 a and multi¬ 
plied this by the number of bins within the ±3 a window. 
All methods gave similar results. The final estimate of 
the number of observed normalizing events was taken as 
the average of the three estimates, (3.050 ± 0.023) x 10®, 
with the uncertainty defined as half the difference be¬ 
tween the largest and smallest estimates (corresponding 
to about Icr, assuming Gaussian statistics). After cor¬ 
recting for acceptance and selection efficiencies, and ac¬ 
counting for the measured branching fractions, the to¬ 
tal number of baryons exiting the collimator was 
(20.2 ±0.8) X 10®. 

Signal-mode branching fractions were obtained from 
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FIG. 5: Observed pK~Tr~ invariant-mass distribution. The 
number of events in the —>■ AK~, A —> pn~ peak was used 
to normalize the upper limits presented here. The arrows 
indicate ±3cr in mass resolution. 

Here N denotes the number of events observed and A 
the acceptance for a given mode, with subscripts sig des¬ 
ignating the signal mode in question and norm the nor¬ 
malizing mode —>■ AK~,A —> p'K~. The number of 
signal events observed is zero in each case. The mea¬ 
sured branching ratios entering into Eqs. (2) and (3) are 
AK-) = (67.8 ± 0.7)%, B{A pn-) = (63.9 ± 
0.5)%, and B{n- E^n-) = (23.6 ± 0.7)% [4]. Sys¬ 
tematic uncertainties include the uncertainties of these 
branching ratios, that of the normalizing-mode back¬ 
ground subtraction, and, most importantly, those associ¬ 
ated with acceptances and event selection. 

To derive 90%-C.L. upper limits for the numbers of 
events observed (17„), we included systematic uncertain¬ 
ties [11] as follows: 

Ur, = Ur,o[l + {Uno-n)a^j2]. (4) 

Here Uno represents the statistical limit based on the 
Poisson distribution with no systematic uncertainties, n 
is the number of observed events, and is the relative 
systematic uncertainty. In our case (n = 0), Eq. (4) 
reduces to 

Uo = 2.3(1-h 2.3fT2/2). (5) 


acceptance, selection, normalization, and branching ra¬ 
tios) was 16% for ^ S°7r“, —> p7r“and 32% 

for —> Att~ , A —> pTT~. Upper limits on the num¬ 

bers of observed events at the 90% C.L. thus determined 
are 2.4 events for > p7r“and 2.6 for —> Att", 

comparable to the 2.3 events obtained from a frequentist 
statistical treatment of the Poisson fluctuation alone [4]. 
We thus obtain B{VL~ —> Att”) < 2.9 x 10“® and 
B{E^ P'!T~) < 8.2 X 10“®, both at 90% C.L. These 
results represent improvements by one to two orders of 
magnitude over previous measurements. 
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